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We investigate the influence of the nucleon-nucleon collision profile (probability of interaction 
as a function of the nucleon-nucleon impact parameter) in the wounded nucleon model and its 
extensions on several observables measured in relativistic heavy-ion collisions. We find that the 
participant eccentricity coefficient, e*, as well as the higher harmonic coefficients, e^, are reduced by 
10-20% for mid-peripheral collisions when the realistic (Gaussian) profile is used, as compared to the 
case with the commonly-used hard-sphere profile. Similarly, the multiplicity fluctuations, treated 
as the function of the number of wounded nucleons in one of the colliding nuclei, are reduced by 
10-20%. This demonstrates that the Glauber Monte Garlo codes should necessarily use the realistic 
nucleon-nucleon collision profile in precision studies of these observables. Our study is relevant for 
hydrodynamic calculations which use the Glauber-model initial conditions. The Gaussian collision 
profile is built-in in GLISSANDO. 
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I. INTRODUCTION 

The wounded nucleon model [l[ and its extensions 
[3, Q have become the basic tool in the analysis of rel- 
ativistic heavy-ion collisions. The simple idea, based 
on the Glauber approach 0, [H| and finding its foun- 
dation in the Landau- Pomeranchuk effect [6[, leads to 
proper description of the bulk properties of the data. In 
particular, the centrality dependence of the multiplici- 
ties from the RHIC energies {^/snn — 42 GeV) [7] to 
LHC {y/SNN = 2.76 TeV) [§] is properly reproduced in 
the mixed model 0], appending the wounded nucleon 
model with some particle production from binary colli- 
sions. Moreover, in the collider experiments it is custom- 
ary to determine the number of participants as a function 
of centrality with the help of the Glauber Monte Carlo 
simulations [1, H-Ul]- 

The Glauber Monte Carlo calculations have also been 
used to determine the initial eccentricity of the system. 
This is a very important quantity, as it gives rise, via 
subsequent dynamical evolution of the system, to the 
observed elliptic flow coefficient, V2- Accuracy of 10% 
or better is requested in this kind of studies to probe 
with sufficient accuracy the sensitivity to various physi- 
cal effects, in particular to the equation of state |13l - [l6| 
or the viscosity parameters [13-26] of the hydrodynamic 
medium. 

The statistical nature of the Glauber Monte Carlo un- 
avoidably leads to fluctuations in the distribution of the 
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wounded nucleons, which yield the event-by-event eccen- 
tricity fluctuations related to the fluctuations of the ellip- 
tic flow [27rt43j , the multiplicity fluctuations as the func- 
tion of the number of wounded nucleons in one of the 
colliding nuclei (the NA49 experiment setup [44]), and 
the overall size fluctuations [45], which can explain the 
magnitude and centrality dependence of the measured 
transverse-momentum fluctuations [46l - l6l| . 

The concept of wounded nucleons can be applied at 
non-zero rapidity, where it has been beautifully tested 
for the deuteron- nucleus collisions [g^], with the conclu- 
sion that the rapidity proflle of the particle emission from 
the wounded nucleons is asymmetric, peaked in the for- 
ward hemisphere, but with substantial emission to the 
backward direction as well. Extension of this approach 
to the nucleus-nucleus case [H, [63] leads to a proper de- 
scription of the forward-backward multiplicity fluctua- 
tions at RHIC (63 - [67j and to a natural explanation of 
the sj_gn and magnitude of the directed flow coeflicient, 
vi 68]. Furthermore, theoretical studies of the correla- 
tions of the forward and backward elliptic flow with fluc- 
tuations and the asymmetric emission profile lead to the 
torque effect [6^, where the angle between the forward 
and backward principal axes of the particle distributions 
in the transverse plane fiuctuates substantially on event- 
by-event basis. All the above-mentioned features show 
that the Monte Carlo implementations of the wounded- 
nucleon approach result in rich predictions and allow for 
numerous dedicated analyses in heavy-ion physics. 

In the wounded nucleon model the basic entity is 
the nucleon-nucleon collision profile, p{b), defined as the 
probability of inelastic nucleon-nucleon (NN) collision at 
the impact parameter b. It is normalized to the total 
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inelastic NN cross section, 



2'Kbdbp{b) — Cinel- 



(1) 



For the most basic quantities, such as the average number 
of wounded nucleons in a nucleus-nucleus collision at a 
given centrality, the shape of the NN collision profile is 
irrelevant. For that reason, in most Glauber Monte Carlo 
codes it is assumed for simplicity that p{b) is simply given 
the step function, p{b) = <d{R—b), with ttR'^ = fiinci- We 
refer to this as to the hard-sphere approximation. 

The purpose of this paper is to show that for cer- 
tain important heavy-ion observables (eccentricity coeffi- 
cients, multiplicity fluctuations, triangular deformation 
[7(]| - l72| ). the shape of the NN collision profile is very 
much important. In particular, the use of a realistic 
p{b), as determined from the low-angle scattering [73| . 
yields noticeably lower eccentricity parameters and mul- 
tiplicity fluctuations, with the reduction at the level of 
10-20% for the mid-peripheral collisions. By providing 
the initial condition, our study is relevant for "preci- 
sion" hydrodynamic calculations, probing the equation 
of state, viscosity parameters, etc. 



II. BASIC ELEMENTS OF THE MODEL 

The unitarity condition links the elastic NN scattering 
amplitude, tei(b), with the inelastic collision profile (we 
neglect, as is customary, the small real part of the NN 
scattering amplitude) 



toi(fo) = 1 - Vl-p(&). 



(2) 



Through passing to the momentum representation one 
obtains 



d^be'^H^lib). 



(3) 



Further, from the optical theorem, the total NN cross 
section is 

atot = 2Tci(0), (4) 
the total elastic cross section is the difference 

O-cl = CTtot - Cincl = j (fb\tci{b)\'^, (5) 

and, finally, the elastic differential cross section is given 
by 



(6) 



with t — —q^ 

In Ref. [Z^ the ISR experimental data [TS - iTSj for the 
total and elastic differential pp cross section were prop- 
erly parametrized with a combination of Gaussians (in 




FIG. 1. (Color online) Comparison of the fit of Ref. ^ 
(dashed line), formula ([7]l (solid line), and the hard-sphere 
profile (dotted line) for the probability of the inelastic NN 
collision, pip), plotted as a function of the NN impact param- 
eter (cTinci = 36 mb). 



that quark-diquark picture it is also possible to repro- 
duce the diffractive minimum in the ISR data). Here we 
use for simplicity a single Gaussian form: 



p{b) = Ae-^^'>^l< 



(7) 



The profiles are compared in Fig. [TJ By construction, all 
fits give the same total inelastic NN cross section. For 
the fit of Ref. ^73^ the elastic NN cross section is 7.8 mb, 
Eq. ([7|) yields 6.8 mb, while the hard-sphere profile gives 
7.4 mb. 

The small difference between the fit of Ref. [t^ (dashed 
line) and formula ([7]) (solid line) is innocuous for the 
study of the heavy-ion observables carried out in this 
work. The parameter A in Eq. (jT]) depends weekly on 
the collision energy, hence we have set 



A = 0.92 



(8) 



for all the studied collisions. The hard-sphere profile, 
shown in Fig. [T] with a dotted line, is, of course, much dif- 
ferent from the Gaussian profile. While the shape of p(&) 
does not affect certain observables, it leads, as shown in 
the following Sections, to noticeable effects in the descrip- 
tion of the shape of the initial fireball or in fluctuations 
of various quantities. 



III. RESULTS OF THE MONTE-CARLO 
SIMULATIONS 

The presented results have been obtained with 
GLISSANDO [ll. The nuclear distributions are generated 
according to appropriate Woods-Saxon densities. The 
correlations from the NN repulsion are generated in a 
standard way by precluding the centers of nucleons to be 
closer to one another than the expulsion distance d. We 
use d = 0.9 fm, which reproduces properly [79| the real- 
istic central correlations implemented in Refs. [sol. Isij. 
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FIG. 2. (Color online) The event-averaged values of the az- 
imuthal eccentricity e (two lower lines) and the participant 
eccentricity e* (two upper lines) for the hard-sphere collision 
profile (solid lines) and the Gaussian collision profile (dashed 
lines), evaluated in the wounded nucleon model for Pb-|-Pb 
collisions and plotted against the total number of wounded 
nucleons, Nw- The panel (a) corresponds to (Tinei = 32 mb 
(top SPS energy) and the panel (b) to Uinei = 63 mb (LHC) . 



A. Eccentricity 

The azimuthal eccentricity is defined as the second 
Fourier moment of the distribution of the wounded nucle- 
ons in the transverse plane. One defines the eccentricity 
coefficient, e, with respect to the axes defined by the re- 
action plane (the coordinate x lies in the reaction plane, 
and y is perpendicular to it), 



Since the principal axes of the distribution fluctuate 
event- by- event, a more appropriate measure is the value 
of the so-called participant eccentricity, evaluated with 
respect to the principal axes in each event [sgJ . Then 



t2 -rr2V 
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xy 



az and a 

y 



xy 



(10) 



are the variances and covariance of 



(9) 



where 

the nucleon distribution in a given event, respectively. 

In the description of the evolution of the system involv- 
ing hydrodynamics [l^ - [26l [82-88]. the initial azimuthal 
eccentricity of the system is a very important quantity, 
as it influences the value of the elliptic flow coefficient, V2 
[s^. Moreover, it is believed it can be used to probe the 
dynamical properties of the system, such as the equation 
of state or the viscosity parameters il7i - t26i] . Thus 

the initial azimuthal eccentricity, supplied as the initial 
condition to the hydro codes, is an important ingredient. 
The fact that e* > £ makes it easier to develop larger 
W2, thus the hydrodynamic evolution may be somewhat 
shorter. This, it turn, is favorable for the uniform de- 
scription of the femtoscopic variables [T^ . 

From the methodological point of view, in order to 
study the details of the system dynamics one should start 
from the most realistic initial condition. Certainly, it is 
not clear what exactly should be used here, as we are 
still far from explaining the earliest-stage dynamics of the 
system formed in relativistic heavy ion collisions. Other 
approaches, such as the Color Glass Condensate (CGC) 
[9(3 . 91] provide a different initial condition [s^ [9^ for 
hydrodynamics than the Glauber models. However, hav- 
ing chosen the Glauber approach, we should use the most 
realistic ingredients. There is no reason to use the sim- 
plistic hard-sphere wounding profile, since the use of the 
realistic profile (as inferred form the pp scattering data) 
is equally simple, while the differences certain for observ- 
ables at the level of 10%. 

In Fig. [5] we show the results of GLISSANDO ]^3] for e 
and £*, obtained in the wounded nucleon model for the 
Gaussian wounding profile of Eq. (O (dashed lines) and 
the hard-sphere profile (solid lines). We note that the use 
of the realistic (Gaussian) profile rather than the hard- 
sphere profile results in lower values of the eccentricity 
parameters. The reduction for intermediate centralities 
is at the level of 10% for the SPS case of iTinci ~ 32 mb 
and at the level of 20% for the LHC case of crinei = 63 mb. 
In fact, the reduction due to the realistic profile brings 
down e* more-less to the value of e evaluated with the 
conventional hard-sphere profile. Thus, the enhancement 
due to the participant geometry is "lost" by the use of 
the proper wounding profile! 

The described quenching effect is easy to understand 
on geometric grounds. Unlike the hard-sphere case, the 
Gaussian profile has a tail extending to large values of 
h (cf. Fig. [Ij. Therefore there is a certain probability 
that a nucleon from nucleus A wounds a nucleon from 
nucleus _B which is "far away" in the transverse plane. 
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FIG. 3. (Color online) Same as Fig.[2]for the case of Au+Au 
collisions with (Tinci = 42 mb (top RHIC energy). Panel (a): 
wounded nucleon model, panel (b): the "hot spot + F model". 



Since there are more nucleons in the a:;-direction (parallel 
to the reaction plane) than in the y-direction (perpen- 
dicular to the reaction plane), with the Gaussian profile 
more nucleons in the x direction are wounded, hence the 
reduction effect. 

Next, we check how the effect depends on the vari- 
ant of the Glauber-like model of the early phase of the 
reaction. The wounded nucleon model is the simplest 
model of this kind. Extensions involve admixing parti- 
cles from the binary collisions (the so-called mixed model 
of Ref. 0,0) nil) or overlaying a distribution of particles 
over the sources (here the "source" means a wounded 
nucleon or a binary collision). The model which leads 
to large fluctuations is the hot-spot model [1, HI]. It 
assumes that the cross section for a semi-hard binary 
collision producing a hot spot is small, (Tbin — 2 mb, 



FIG. 4. (Color online) The event-by-event scaled standard de- 
viation A£*/e* for the Au-j-Au collisions with (Jinci = 42 mb 
(top RHIC energy), evaluated for the hard-sphere profile 
(solid line) and the Gaussian profile (dashed line), plotted as a 
function of the total number of wounded nucleons, Nw- The 
horizontal line shows the limit of 1^4/71 — 1 derived in [3^ . 



but when such a rare collision occurs it produces a large 
amount of transverse energy equal to atTw/crbin- Then, 
we overlay the F distribution on the produced sources, as 
described in Ref. Q. The purpose of using this model is, 
as mentioned, to investigate a case with large statistical 
fluctuations. 

The results are shown in Fig. [3l where we show the ef- 
fect of the reduction of the aziniuthal eccentricity for the 
case of the highest RHIC energy, comparing the wounded 
nucleon model (top) and the hot-spot-j-F model (bot- 
tom). We note that for both cases the reduction is sim- 
ilar, hence the effect does not depend on the particular 
Glauber-like model used for the analysis. 

In Fig.|4]we show a measure of the event-by-event fluc- 
tuations of e* , namely, the scaled standard deviation of 
this quantity, investigated in Refs. [27i - i33l . Issi , [37l - l4l| in 
connection with the fluctuations of the elliptic flow. We 
note a small increase of the fluctuations for this quantity 
when the Gaussian collision profile is used compared to 
the hard-sphere case. 



B. Higher azimuthal Fourier components 

A similar reduction effect occurs for the higher Fourier 
components of the azimuthal distribution of sources. The 
coefficients of rank defined in each event as 

[3i,[73 

, E.Pfcos[n(^.-0*)] 
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tation of various observed phenomena, such as the ridge 

mill- 
In Fig. Owe show the wounded-nucleon predictions for 
e* with 71 = 3, 4, and 5. Similarly to the n = 2 case, we 
observe a sizable reduction when the Gaussian collision 
profile is used. For mid-peripheral collisions it is at the 
level of 15%. Therefore the realistic NN collision profile 
should be included in the analyses of the triangular and 
higher-harmonic flow using the Glauber approach. 



C. Multiplicity fluctuations 

Another important class of phenomena studied in rel- 
ativistic heavy-ion collision are the multiplicity fluctua- 
tions [13, [13, [13, which carry information on the NN 
correlations. We have tested several reactions, in par- 
ticular those to be investigated in the NA61 experiment 
[99.- 103]. The basic measure of the multiplicity fluctua- 
tions is the scaled variance of the wounded nucleons in 
the target. 



var(Ntarg) 
(^targ) 



(13) 



plotted as a function of the wounded nucleons of the pro- 
jectile, which are measured experimentally in the NA49 
or NA61 setup. 

A typical result is shown in Fig. [5] for the Ar-fCa col- 
lision with (Tinci = 32 mb. We note a reduction of Wtarg 
when the Gaussian profile is used compared to the hard- 
sphere case. The effect is at the level of 15%. A similar 
effect is found for other combinations of the colliding nu- 
clei. 



D. Other results 



FIG. 5. (Color online) The wounded-nucleon model predic- 
tions for the third ( panel (a)), fourth (panel (b)), and fifth 
(panel (c)) Fourier coefficients of the azimuthal distribution. 
See the text for details. 



where i runs over the number of sources (wounded nucle- 
ons or binary collisions) in each event, pi is the distance 
of the source from the center of mass of the fireball, and 
is its azimuthal angle. The angle (f>* is adjusted in 
each event in such a way as to maximize e* , which gives 
the condition 



tanfli* = 



E» Pl cos(n0i) ' 



(12) 



For n — 2 the above definition overlaps with Eq. ([TT 
The odd rank coefficients arise entirely due to statistical 
fluctuations [3^ . 

The hydrodynamic evolution, when started from initial 
distributions carrying e* , generates the collective flow co- 
efficients Vn (tiI . [7^ . They are important in the interpre- 



We have also carried out other analyses done in heavy- 
ion physics, using the Gaussian profile in GLISSANDOjj. 
The forward-backward multiplicity correlations 0, Il04| - 
llQ6l | are not affected by the choice of the NN collision 
profile. The same is true of the size fluctuations studied 
in Ref. (45j . which influence the transverse- momentum 
fluctuations. The size variable is deflned in each event as 



51 A' 



(14) 



where i runs over the sources and pi is the distance from 
the center of mass of the fireball. The result for the 
event-by-event scaled standard deviation of r is shown in 
Fig. [T] We note a very small effect from the choice of 
the NN collision profile, with a small reduction when the 
Gaussian profile is used. 



IV. CONCLUSIONS 



We conclude be reiterating our main results: 
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The Glauber Monte Carlo codes for relativistic 
heavy-ion collisions should definitely use the realis- 
tic NN collision profile, as it affects important ob- 
servables in a sizable way. 



The eccentricity parameters, e*, and the higher 
rank parameters, e* , are reduced by 10-20% when 
the realistic (Gaussian) NN collision profile is used, 
compared to the case with the hard-sphere profile. 
The multiplicity fiuctuations, as measured in the 
NA49 or NA61 experimental setup, are reduced by 
10-20% when the realistic NN collision profile is 
used compared to the hard-sphere case. 



FIG. 6. (Color online) The multiplicity-fluctuations measure, 
ijJts.is, for At+Csl collisions with ctiuci = 32 mb, evaluated 
for the hard-sphere collision profile (solid line) and the Gaus- 
sian collision profile (dashed line), plotted as a function of 
wounded nucleons in the projectile, N^"' . 



• The initial condition following from the realistic 
wounding profile should be used in hydrodynamic 
studies. 
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We note that GLISSANDO [1(3 has a built-in realistic 
(Gaussian) NN collision profile and can be readily used 
for those analyses where sensitivity to the NN collision 
profiles appears. 



One if us (WB) wishes to thank Adam Bzdak for useful 
discussions. 



FIG. 7. (Color online) The event-by-event scaled standard 
deviation of the size variable r for Au-f Au collisions at (Jinoi ~ 
42 mb (top RHIC energy) for the hard-sphere collision profile 
(solid line) and the Gaussian collision profile (dashed line). 
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